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SUMMARY 


■A method is given for calculating the temperature of a 
surface wetted either by a pure liquid, such as water, or by 
a mixture, such as alcohol and water# The method is applied 
to, the problem of protecting, by alcohol, propellers and the 
induction system of the engine against ice, The minimum 
quantity of alcohol required is calculated for a number of 
arbitrarily chosen conditions. The effect of evaporation of 
alcohol is shown by repeating the calculations for a nonvola' 
tile fluid. The method can be applied to other problems in 
evaporation, for instance, to the evaporation of fuel in the 
induction system of the engine. The psychromstric equation, 
used in' wet-*bulb hygrometry, is deduced in its general form. 
The effect of kinetic heating is included in this equation. 


INTRODUCTION 


The subject Of this note is the evaporation from a wet 
surface in a stream of air, the surface being unheated other 
than by convection of heat from the air. A familiar example 


♦This report was prepared by Mr, Hardy in collaboration 
with the staff of the Ames Aeronautical Laboratory during a 
period of active participation by Mr, Hardy in the NACA icing 
research program-. 
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is the wet— hulb thermometer, •iTapdration is considered pri- 
marily in connection with the use of volatile organic liquids, 
such as alcohol, to prevent the formation of ice on parts of 
an aircraft, such as the propeller, or the induction system 
of the engine, A method is presented by which the temperature 
of a surface, which is wetted' by a liquid consisting of one 
or more volatile components, may be calculated. It is pos- 
sible, therefore, if the rate at which water reaches the sur- 
face is known, to calculate the minimum rate at which any par- 
ticular de-icing fluid must be supplied in order to prevent 
freezing. ■ r,:;,-'. 

In the case of a surface wetted with water, the equation 
derived is the psy.chrometrio equation in its complete form. 
Since the effect of kinetic heating is inciudeds this equa- 
tion may be used In the measurement of humidity, by' the vret— 
bulb method, at high speeds. This equation is of use, also, 
in calculating the rate of heating required to protect a 
surface exposed to conditions of icing. The datum tempera- 
ture, in calculating the rate ©f dissipation of heat from 
the surface by the method given in references 1 and 2, is 
the wet— kinetic temperature of the surface. This may be 
calculated from the psychr bmetric equation, in the form giv- 
en in this report, as ah alternative to the method proposed 
in reference 3, The use of the specific heat of wet air, as 
proposed in this- reference, is open to the objection that, it 
appears to require that the saturation of the air, in the 
vicinity of the surface, is maintained by evaporation of’ the 
droplets of water carried by the air. It appears, however, 
from the different approach to the problem presented in this 
report, that it is unnecessary to postulate a particular 
mechanism of evaporation. 

This report is complementary to that of reference 1, 
since it deals primarily with the theoretical aspects of 
protection against ice by chemicals. The methods of cal- 
culation, though 'applied to the problem of icing, are appli— 
cable to problems of evaporation generally. They could be 
used, for Instance, to calculate the rate of evaporation of 
■fuel from a heated surface in the- induction system of the 
engine, 

SYMBOLS 

The following symbols are used throughout this report'; 
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N rate of evaporation per unit area, pounds per second, 
sq.uare foot 

H rato of transfer of hoat per unit area, British thernal 
unit per second, s'q..uare foot 

iCg coefficient of evaporation, dimensionless 

coefficient of evaporation of water, dimensionless 

kj^ coefficient of transfer of heat, dimensionless 

coefficient .of surface friction, dimensionless 

J mechanical equivalent of heat, foot-pounds per British 
thermal unit 

g gravitational constant, feet per second squared 
Vq velocity of free stream, feet per second 

V 3 _ velocity locally at edge of houndary layer, feet per 
second 

p density of air, pounds per cubic foot 

pg density of vapor, pounds per cubic foot 

M molecular weight of air (- 28,96) 
a 

molecular weight of water vapor ~ 0,622) 

Mg molecular weight of vapor 

n concentration of vapor, pounds per pound of air 

e vapor pressure, millimeters of mercury 

e' vapqr pressure saturated air at temperature t’, 
millimeters of mercury 

p Barometric pressure of air, millimeters of mercury 

(j, viscosity of air, pounds per second, foot 

k thermal conductivity of air, British thermal units per 
second, square foot degree Fahrenheit per foot 
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specific heat of ,air , Br it Ish ■ therraa/l -anlts per pound, 
degree Pahrenheit 

specific heg,t of wet air, British thermal units per 
" pound, degree Pahrenhelt 

D diffusivity of vapor in air, foot squared per second 

L latent heat of vaporization, British thermal units per 

pound 

t temperature , degrees Pahrenheit (unless otherwise 
specified) 

t' wet— hulh temperature , degrees Pahrenheit- (unless other- 

wise specified) 

Tr Taylor’s number, dimensionless. , • 

Pr Prandtl 's number, dimensionless 


Subscripts 

0 refers to ambient air at static temperature 

1 refers to conditions at the edge of the boundary layer 

s refers to the surface 

w refers to water vapor 

TH20EY 

In the analysis which follows, it is assumed that the 
surface wetted by the f luid is isolated t her raa lly , so- that 
it does not either receive or lose heat by conduction or 
radiation, I,n these, circumstances , the. temperature assumed 
by the surface is such that the heat lost by; eva.po.ratlon' 
equals that gained by convection. 

Rate of Svaporation, 

The rate of evaporation f.rom unit area of a surface in 
air flowing at velocity 7 is given by the equation 
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( 1 ) 


In this is the coefficient of STapor at ion w hich is non- 

dimensional ^ and which is related by EeynoXds* analogy with 
the coefficient of surface friction The potential in 

the process of transfer is ^ ^ difference, namely, 

between the concentration of vapor in equilibrium with the 
liquid at the surface, and that in the air outside the bound- 
ary lawyer « The concentr at ion is expressed as the mass of 
vapor per unit mass of air 5 it may be expressed in ’terms of * 
vapor pressures by substitution from the equation 

M 0^ 0 

n 

“a 

Equation ( l) then he comes 



( 2 ) 


in which e^ is the Tapor pressure in equilihrium with the 

liquid at the temperature of the surface the saturation 
vapor pressure , that is, in the case of a pure liquid. 

Equations (l) and (3) can he applied irrespective of 
whether the process is one of evaporation or condensation. 

The value of k^, the coefficient of evaporation, is 

determined, fundamentally, hy the thickness of the houndary 
layer and the nature of the flow, as defining the field of 
diffusion, and also hy Taylor's number, a parameter which 
gives a measure of the activity of the process of diffusion, 
Taylor's number, Tr , is defined as 


Tr 


M. 

Dp 


in which D is the diffusivity of the vapor, and y, and ' p 
the viscosity and density of the gas through v/hich the vapor 
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is diffusing. Taylor '' s' ntimler, .is. the r'atio.,, of thfe, d'iff u.s iv--. 
ity" of jn'omentuin , usually called hinemat'ic viscosity, to the 
diffusivlty of the vapor. ' The eiuat-io''h' whle'h'defin'es' d'lf--'''' 
f usivity is • 

1 dm ^ ' - ' V e 
,A dl -dy 

in which dm/ dT is the mass of vapor diffusing in unit time 
across an area A, and d^, / dy is the concentration gra— 

dieht . When applied to the , d iffus i on of momentum, in the 
gas through which the vapor ' is ' diffusing , this eq.uation he — 
c 0 me s 

_ i = tpv* Me.Il ' ' ' ■ ' ■ ' 

A dl ~ " dy 

In this eq.'uation (P ) ' • is the " s.e If — dif f 'us ivi ty" of the gas, 
and from the kinetic theory of gases this is equal to d/p • 

On suhstituting for (D)’, the equation hecomes the normal 
equation: .whi'ch defines .v.iscosity. Taylor ' S ■ mitoher. isths'.’ .. 
rat ip , of . _ (Dl 'to . .• ,,,t. , 

The conduction of heat, also, i's a process of diffusion — 
of the energy, of -molecular motion.; The dif fus iyity - of .heat 
in air-; is,,, k/p Cp and,, with . appropriate : suhst i tut ions., equa- 

tion (i) hecomes the normal equation which defines the conduc- 
tivity of' heat.. . In the case o.f.heat, the ratio of. -the; dlffu- 
sivity of momentum to that of heat is called Prandtl's miraher. 



The value ' of Tayl or ’ s ' nuither like Prandtl ' s ■ numher", is - 
independent of pres sure , hut cha'ngQS' s lowly ' wi t'h temperature-. 
T-he- dif f usivity of a vapor in te-rms of its ■value" Dg.j. , at ' ' 

0° 0 and a pressure of 1 atmospheie is given* hy the equation 



Ps t 

'^p'" 



wfiere the temperatures are in degrees ahsolute , ' and ffl ' has' 
a value which' range s'' f f om" 1 . 75 t o .P . 0 . ' ■ 

The coefficient of evaporation k^ is related, thro'agh 
Reynold's analogy, to the coefficient of transfer of hea-t 
as defined hy the equation 
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H 



Atj_) 


( 5) 


In this equation, At^^. is. the iucrement in t emp eratur e , fr am 
kinetic heati.ng, which is discussed-in a later section. . 

Reynolds analogy is, hetwedn heat ' transf er and surf ace 
f r ict ion 5 ^ arid equht ions. . r elat ing- the 'coefficients are de— 
rived'in re.fer.ence .4, pages 623 to. 626 , and 654' to 657,. 
These equat ions , with suhs titut ion of coefficients and of 
Tr for Pr, apply also to mass transfer. For instance, 
when the flow is laminar 



When the flow is turhulents the equation corresponding to 
equation (6) is 

(7) 

This equation is approximate,- hut is- sufficiently accurate 
v;hen Pr ' has a value near to 1, A more accurate equation, 
given hy Von K^rmdn, is 

i + 5 /X - 1) + log fl + I (Pr - l)l 

Cf J ^ \ ® jj 

The preceding equations, although derived for the flow across 
a flat plate and through a pipe, can he applied quite gener- 
ally and without restriction to determine the value of kg 
from that of kj^. 


Jn. 

^^h 


+ 5.6 ( Pr - 1) 


(2 

C. 
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Teaperature of a Wet • Surface 

The temperature assumed by a \iret surface which is thor~ 
mally isolated is such that the loss of heat by evaporation 
equals' the gain of heat by convection from the air. The 
surface is isolated thermally in the sense that there is no 
gain or loss of heat by conduction or radiation.' For such 
a s'urfa,ce the balance of flow of heat is given by the eq^'da'^- , 
t ion , ' ' 

: , L¥ = H - , . ( 8) 

One volat ile c6mp one nt For a surface wetted either by 
a pure liquid or by a solution in which one component, only, 
is volatile, substitution in equation (8) from equations (3) 
and ( 5 ) gives , 


I* ^e ^0 P 


Me 

“a 




Pi- e- 


J 


- kh p (t3 


At ) 


This may be rearranged for convenience in calculat ing th e 
temperature of the surface as follows: 


t 


ke Me 

t^ + Af = -- ~ 

K M 

h a 


(: 



( 9 ) 


An approximate form of this equation, which may be used at 
low temperatures when the vapor pressure is small relative 
to the barometric pressure, is 


- ,^s + ^^1 = 



M e 
e /^s 


M 


a 


(■ 



) 



( 10 ) 


T w o _ V 0 1 ai j, c n e , — ¥ h g n t h e s u r f ace i s e 1 1 e d hj 

a liquid consisting of Tolatile components, suck as al- 

cohol and watpyrk each evaporates indepondont ly of the other, 
so that the effect is additive • For a system of two comp^o — 
nonts, equation (10) becomGS 
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(li) 


The second component is indicated hy the change in suhscript. 
The vapor pressure of each component in eq.uilihrium v;ith the 
surface.s , Og a.nd that appropriate to the mixture. 

It differs from the vapor pressure of each component in the 
pure state. It. is important to notice that it is the compo- 
sition of the liq.uid at the surface which is important. 

This, in general, will differ from the composition of the 
hulk of the liciuid unless there is mechanical mixing, since 
the rates of evaporation of the components will he uneq.ual; 
in fact, under certain conditions, one may he evaporating 
while the other is condensing on the surface. 


Part ia ll y w e tted s u rfac e.— Equations (9) and (lO) ap”^ly 
to a surface which is completely wetted with liquid. Condi- 
tions may he such that less liquid is received hy the sur- 
face tha,n coxild he evaporated. This condition of .partial 
i^etness requires that the equations he modified hy multiply- 
ing the right-hand side hy the wetness fraction. The wetness 
fraction is determined hy dividing the rate at which water 
is received hy the surface hy the rate at which it v;ould he 
lost hy evaporation if the surface were complotoly wetted. 

The rate of eve.poration is calculated from equation (.3). 


The psych r ometr ic eQuat ion.— The psy chrome tr ic equa.tion, 
quoted later as equation (18), is the empirical expression 
of a considerahle body of experimental data. It applies only 
at rates of air flov.', past the hulhs of the thermometers, so 
that kinetic heating is ina.ppr e c iah 1 e , and so that the field 
of pressure around the wet hulh is of negligihle intensity 
as compared with the harornetric pressure in the free stream. 
At high rates of flow, it is evident that the shape of the 
hulh relative to the direction of flov/ of air is of grea.t 
importance. For instance, in the case of a cylinder across 
wind the temperature tg as given hy equation (9) will vary 
around the cylinder because of variation of At, and of p ^ 
and e^. If the cylinder is of material of high thermal 
conductivity, the temperature assumed will he such as to 
satisfy the balance of flov/ of heat over the whole surface 
of the cylinder. The temperature of the cylinder, which v/ill 
ho the apparent wet hulh temperature, will ho such as to sa.t— 
isfy the following equation 


•'h P '^^0 *^p ^ ‘'h 


-ti-Ati)ds = Z' k^. 


M 


V w/®h 
P ^0 M~ \ “ 

'•Pi - e 1 


^w 'I s 
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in which t-j^ is the temperature of the cylinder, and e-|j 
the vapor pr e's suf’o ^ f o^ saturation at this temperature. 

In psychr ome t r y.,’ what is req.uired is the vapor pressure 
of the undisturbed stream, e^, and it is evident that this 
cannot he obtained- hy any simple process of calculation from 
the foregoing equation. It can be obtained, however, if the 
bulb of the thermometer is in the form of a flat plate, edge 
on to the flow of air, or is such as to satisfy the require- 
ments of cons.tancy both of kinetic temperatu.re and of pres- 
sure over the, surface. 

In the, case of a thin flat plate, both barometric and 
vapor . press ur'es at the edge of the boundary layer are iden- 
tical with those in the free stream. Equation (9), there- 
fore, may be rewritten as 


tg + At^ = 


k,, M,- 
w w 


^h \ 


c?, 


^0 

*0 ^ 0 ' 




( 12 ) 


The value of '^31 this equation, is the temperature of 

the wet plate, and the equation may be used for psychrometry 
at high speeds. 

Kinetic Heating 

The effect of kinetic heating has been introduced into 
the equations, which give the balance of flow of heat to and 
from the surface, by the inclusion of At^, the increase in 
temperature caused by kinetic heating. The temperature of 
the surface when the rate of convection of heat is zero, 
therefore, is t^ + At^. 

The value of At ^ when the flow is laminar is given, 
in reference 4, by the equation 


A t 1 


2g JCp 




( 13 ) 


When the flow is turbulent , reference 5 gives 
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These Gq.uations presume that there is no evaporation 
from droplets of water \i;hich, under certain circumstances, 
may he present in the boundary layer, since the specific 
heat of dry air is used. This is consistent w'ith the assump- 
tion, in equation (8), that evaporation occurs only at the 
surface. The effect of evaporation within the boundary layer 
will he discussed later. 


Outside the boundary layer, the effect of a change of 
phase from vapor to liquid, is of direct importance, as it 
affects the value of t^, and also the value of e^. As air 
flow round a body, and its velGcity changes from to 7^, 

the process is isentropic. Tho value of ti, on the assump- 
tion of no change of phase, is. 



and tho value of is 

e>i 



(15) 


Over the ferward part of an airfoil, or other body, it is 
probable that there is no change in phase even in air which, 
initially, is saturated. At the other extreme, it may be 
assumed, in the case of an increase in velocity, that con- 
densation occurs so rapidly that tho air remains saturated. 
In this case, the value of t^ is given by substituting the 
specific heat of wet air, equation (15), the value 

of Cp^ being givexi by equation (17). The value of e^ is 
the vapor pressure for saturation at t^. 3o^us-tion (17) can 
be modified to suit a condition of partial condensation >;ith 
super saturat ion, a condition intermediate between the ox- 
tromos just defined. 


Kinetic tempera ture in cloud.- When droplets of water 
are carried by the air stream, as is the case v/hen an ,.air — 
craft passes through cloiid, the method proposed in refer- 
ence 3 for calcula.ting kinetiG temperature is to use tho 
spocific heat of wet air in equations (13) and (14). The 
objection to this is that it appears to require evaporation 
from the droplets sufficient to- maintain a co.ndition of 
saturation throughout the boundary layer. It can be shown, 
hovrever , -that it is the evaporation at the surface v/hich is 
of primary importance and, in th-e case of \-/ater, that the 
role of the droplets , other, ths-n in wetting the surface , is 
entirely neutral. 
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In the case of water, the diffusirity of the vapor is 
eq.ual to the diffusivity of heat ,■ as, v;lll he shov;n in the 
next ‘section. The distr ihutlo'n of vapor in the boundary 
layer, .therefore, is identical with the distr ib.ut ion of teEi— 
perature ,if each is expressed in the nondiraensional form 

e ™ e ^ t 

e g — s 1 . ■fc g “ 't 1 

If the air' outside the boundary layer is saturated with vapor 
and • the surface is wet, the relation between the distribution 
of temperature and vapor pressure is such that, a condition of 
saturation obtains throughout the boundary layer., It would 
be anticipated, therefore, that the kinetic temperature of a 
surface in cloud, as calculated from the specific heat of wet 
air, would be the same as that calculated from equation (10) 
for a wot surface in clear saturated air. The identity may 
be shown formally. 


The temperature of the surface, when calculated from the 
specific heat of wet air, is 


! ■ t 1 .+ A t 1 


w 


^ 16 ) 


The value of At is given by the right-hand side of equa- 
tion (13) or (14) 'with Cp^ substituted for Op, so that 


At 


Iw 


= At 


Pw 


in which 


My Ly 8 e 


° Pw ° P Ma P a 5 t 


(17) 


The value of 8e/6t is the mean for the interval of tempera— 
t ur 0 , t g — 1 1 , so t ha t 

5e _ eg ~ ® 1 
8 t t g ■*-- t ^ 

By substitution in eq^iation (16) 



1 + 
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When irml tipi led out, this gives 


t J.+ 


At. 


w 


a 




p 


y c 


p 



This differs from eq.uatioji (lO) only in the absence of the 
ratio When this ratio is unity, which ToQ,vAroB 

that the dif f us ivlt ies of heat and of vapor are eq.ua! , the 
eq.uations are identical* 


IDvaluation of Ta.ylor*s Humber for Walter Taper 

There is some uncertainty as to the correct value of 
Taylor *s number for water vapor, owing to uncertainty in 
the value of the diffuslvity of the vapor. In the caseof 
water vapor diffusing through.air, the values for’ diffusiv— 
ity obtained by different exper iinent er s range from 0*198 
to 0.252 centimeter^ per second. The value may be deduced 
from the psychr ometr ic eouation. This is given in reference 
6 a s ^ 

e = e' - 0.000652 p (t - t ' ) (l + 0.00102 t ' ) (is) 


for pressures in millimeters and temperatures in degrees 
centigrade. It should he compared with the accurate form 
given in equation (12). Transposing equation (12), and 
using the same sjmhols as in (18) , with p in pla-ce of 
— "e > gives 


e 


1! C h V 

p _a la. ail (i 

M L k 

V/ W Ui,' 



so that for an air temperature of 32'^ R 



0.000652 




1.004 


It appears, therefore, that Ta.ylor's nuraoor for w^ater vapor 
is almost exactly equal to Prandtl's numher at this temper- 
ature. The value for Prandtl’s numher at 32® P is 0,71. 
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To make Taylor's number eq.ual to this, the dif f tvs ivity of 

water vapor must be 2.01 x 10“^ feet^ per second (0.136 

cm^/sec). The value given in the International Critical 

Tables is 0.22 centimeter^ per second. It is believed that 

the value deduced from the psychr ometr ic equation is the 

more reliable, so that in evaluating equations (9) and (lO), 

for a water— wet surface, the ratio k ;k, should be taken 

w h 

as unity. The ratio changes slowtly with temperature; for 
instance, at 0° J it is 0.996 and at 60*^ F, 1,007, 


RSSULTS OF CALCULATIOHS 

A number of calculations have been made in order to 
demonstrate the application of the preceding equations to 
specific problems in connection with protection of s,ircraft 
against the formation of ice. The kinetic temperature of 
a water— wotted surface has been calculated in order to 
demonstra.to the effect of barometric pressure and tempera- 
ture of tho air. Tho offeset of evaporation, when a vola- 
tile organic liquid is used to prevent the formation of ice, 
is demonstrated by calculations of the minimum amount of 
alcohol required for particular conditions, and comparing 
this with the amount required vtfhon a nonvolatile liquid, of 
othorviriso similar properties, is used. The liquid chosen 
is ethyl alcohol since there are complete data for this in 
the Internationa/l Critical Tables. A calculation has been 
made 8,1s o for methyl alcohol to domonstrate the effect of 
a change in the nature of tho liquid. 

These calculations have been made for a section of a 
blade of a propeller. Calculations have been made also for 
the carburetor to show the mechanism whereby ice forms when 
the temperature of the ambient air is above freezing and 
the quantity of alcohol necessary to prevent the ice. 

In these calculations the data for tho va,por pressure 
of alcohol, both pure and in solution with water, have been 
taken from the International Critical Tables. The value of 
Prandtl 's number has been taken as 0.71 and the value of 
Taylor's number for wa,ter vapor has been taken as equal to 
this . 


Kinetic Tomporaturos of Propeller Blade in Wet Air 
Calculations have been made of the temperatures at a 
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point on the surface of a propeller h la-do , taking the local 
velocity as 600 feet per second a,nd the flow as laminar. 
These have heen made for harometric pi’essures of 760 milli- 
meters and '350 millimeters, and for air temperatures (static 
of 0° R and 25° F, the air being a, ssumed to he saturated in 
each case. The results are given in ta,hle I. As an .in- 
stance' of the method, the, calculation for a temperature of 
0° F and pres. sure of 760 millimeters may bo cited. For 600" 
feet per second, the value of At, for dry air, from eq.ua~ ■ 
tion (13) is 25.3° F, so that equation (10) for' the cond.i— ■ 
tion specified is . . 

0° - t., + 25.3° = (og - 1.11) 

. . , ^'.eo ® 

By trial it is found that a value of 19.4° T for t^ sat- 
isfies this eq.uo,tion. The substantial decrease in kinetic 
temperature with decrease in pressure is of practical sig- 
nificance in the translation of the results of wind-tunnel 
tests, on the formation of 403, mado at ground level, to 
the condition which will occur at high altitude when the 
protection from kinetic heating is less. 

The kinetic temperature of a metal blade of a propGller 
is influenced by the conduction of heat, both along the 
blade and through the blade from the pressure face to the 
cambered face. If tho blade is hollox";, or is of material 
of low thermal conductivity, the lowest temperature will be 
at the position where the local valocity is grea.tost. 

The effect of local velocity is shown by c.alcuiat ions 
for a velocity 1.4 times that of the blade of a. propeller 
v/hieh is taken a,s 600 feet per second; the local velocity, 
therefore, is 842 feet per second. The incroaso in veloc- 
ity choson is -less than occurs on the cambcrt3d face v/hon 
the propeller is working at high incidence. The static 
tomporature of the undisturbed air is taken as 25° F, the 
barometric pressure as 350 millimeters, and tho air as be- 
ing saturated with wa/ter vapor. At the position when the 
velocity roaches 842 feet per second, the pressure, calcu- 
lated by Bernoulli's equation for compressible flov/, is 282 
millimeters. The static temperature of the air, locally, is 
given by equation (15) if it is assumed that the decrease in 
pressure is so rapid that.no condensation occurs in the air. 
A calculation has been made on this assumption (a), and an 
independent calculation has been made assuming complete con- 
densation (b). 
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(a) Assuming no condensat ion: ■ The reduction in temper- 
ature of the air is 28.7° P so that the static temperature 
locally, "by eQ.uation (15), is —3.7° P. The vapor pressure 

is 3 . 43 X = 2.76 millimeters. If it is assumed that 

350 

the flov.’- is laminar , At ^ has a value -of 49.6° P and equa- 
tion (10) gives the value of t^ as 30.5° P. The effect of 

kinetic heating, therefore, is an increase in tomperature, 
locally, of 5'. 5° P; at the stagnation point in clear air the 
increase is 30.0° P. 

(h ) Assuming complete condensat ion: 'The red\iction in 

temperature of the air is given hy equation (16) with the 

specific heat of wet air C_ . suhstitutod for G_ . The 

Pw P 

value of 0 is given hy equation (17) v/ith substitution 
•9V7 ' . ' 

from the equation 

5 e s Q “ e 1 

5t t^ - ti 

in which the subscript o denotes the initia.l valtie , and 
1 , that after expansion to 2 82 mi 1 1 ime t er s pressure. The 
values of t ^ and must be found by trial to be con- 

sistent, approximately, with the value ti calculated from 
equation (15), The reduction in temperature is found to bo 
13.6° P, so that tho static temperat-are of the air locally 
is 11.4° P. The value of tg, which can be calculated 
either from equation (16) or from equation (lO), is 32.9° P. 

In tho foregoing calculations it is assumed that tho 
surface is wot w^ith water oven when the terapcrat^lro is be- 
low freezing. The temporature of the surface so calculatod 
gives the datum temperature from which oither the rate of 
heating to prevent ice, or the. rate of icing, may be calcu- 
lated by the methods of reference 2. 


Blade Temperatures with Alcohel- Water Mixtures 

Calculations have been made for the ease in v/hich ice 
is prevented by supplying ethyl alcohol to the blades of 
tho propeller. The calculations have been repeated for a 
nonvolatile fluid having the same deprosssint effect on the 
freezing point as ethyl alcohol. The difference in tho 
quantity of fluid required gives the excess alcohol which 
must be supplied in order to neutralize the r of r igerat ing 
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effect caiAsed "by Its ovaporat ion . • A similar calculation 
has ‘been made for methyl alcohol to shovr how volatile fluids 
of different composition may he comparod. 

In those calculations the c encentrat ion of alcohol in 
solution ill water is such that' the froozing point is de- 
pressed exactly to the temperature assumed hy the hlade. 

This may he termed the correct mixture for the particular 
conditions. A few ealculatioas have been xaade for a sur- 
face I'/etted hy pure alcohol to illustrate the effect of 
ma,ldistr ihution of fluid .in cooling the hlade hy conduc- 
tion of heat from the warmer to -the colder parts. The 
calculations were worked with the temperature in degrees 
centigrade hecause the data in the International Critical 
Tahles are given in this scale. In many cases the calcu- 
lations have hoon made in reverse, starting vxith an assumed 
temperature of the hlade and calculating the temperature 
of the air roq.uired to produce this. This procedure allows 
a more direct comparison since the effect of evaporation is 
determined primarily hy the temperature of the- surface. 


Calculations have heen raade for a speed of 450 feet 
per second assuming laminar flow, for hlade temperatures of 
2S° R (-50 C) and also 10 . 60 R (-12° C), and for harometric 
pressures of 760 millimeters a,nd 350 millimeters. Two cal- 
culations have heen made for a, speed of 250 feet per second. 
The results are given in table II. As an instance of the 
method, the calculations for the hlade . temporature of 23° R 
at 760 millimeters pressure, the first in table II, may he 
cited. Por this condition, the concentration of alcohol 
reauired to dopress the freezing point to 23° P, from the 
data of the International Critical Tahlos , is 124 graras per 
liter of wa,ter . The v.apor pressure of alcohol in eq.uilih— 
rium with this mixture is 0.172 times the va,lu3 for pure 
alcohol at the same .temperature — na,mely , 1.44 millimeter s 
at 23° P. The vapor pressure of water is 0.95 times the 
value for pure water — namoly, 3.00 millimeters. S-ahstitu— 
tion in equation (11 ) for a speed of 450 feet per second 
gives 




23 + 14.2 = if (1.44 
760 


0 ) + 


(3.0 

760 


e 


) 


The effect of eva-peraticn is found by calculating the 
temperature that the blade v/oxild assume in air ,at 15.2° P 
v/ith a nonvolatile fluid, and calculating the concentration 
of fluid required. The fluid is assumed to liave the saxae 
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depressant effect on the freezing point as ethyl alcohol. 
Substitution in e<luation (lO), since in this case there is 
one volatile component, gives 

15.2 - t„ + 14.2 := (0.96 e„ - 2.23) 

® 760 

The vapon pressure of water in ec^uilihr ium with the mixture 
is 0.96 times that of pure water. From this eq.uation, the 
value of t„ is 25.5° F, and the concentration of fltiid of 

the same raolecular weight as ethyl alcohol is 92 grains per. 
liter of water. The excess of ethyl alcohol roouired to 
noutralizo the refrlgorating effect produced by its evap- 
oration, therefore, is 

100 = 35 nercont 

\ 92 / 

The calculations show that the offoct of evaporation 
decreases with tomporature and increases 'with altitude.. 

In the range of temperature in which ice is usually encoun- 
tered, 15° to 32° F, the wastage of ethyl alcohol on ac- 
count of its volatility is considerable. The nonvolsi. t ile 
fluid with which the alcohol is compared is hypothetical, 
but the concentrations of an actual fluid of low volatility 
may bo calculated from those for the nonvolatile fluid, 
which are presented in table II. A fluid such as ethylene 
glycol, for instance, is practically nonvolatile, and, 
since it does not dissociate in solution, the amount re- 
quired will bo in the ratio of its molecular v/eight to that 
of tho nonvolatile fluid which is the same as that of ethyl 
alcohol. The concentration of the nonvolatile fluid in 
table II, therefore, should be multiplied by 1.35 to find 
the equivalent concentration of ethylene glycol . 

The offoct of speed is shown hy comparing tho concen- 
tration of alcohol required at 450 feet per second, Ko . 2 
in the table, with that required at 250 feet per second 
at tho same temperature of the blade. Ho. 3, and also ’with 
that required at the same temperature of the air. Ho. 4. 

The last Shows the requirements at different stations along 
the blade of a propeller exposed to the same conditions of 
icing. A comparison on the basis of the concentration of 
alcohol may be misleading, because it takes no account of 
the rate at which water is caught by the blade . The rate 
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of catch, roughlj", is proportioual to tho speed, so that 
more alcohol, actually, is required at the outer station 
than at the inner for the ps.rticular conditions for which 
the calculation wa,s made. 

A direct comparison between methyl and ethyl alcohol 
is given by the calculations for an a,ir temperature of 
19.2° E, Hos. 2 and 9 in the table. The weight of methyl 
alcohol, per 1000 grams of w'ater, is slightly loss than 
that of ethyl alcohol. In practice t-he 8,dvantage is more 
than nullified, probably, by maldistribution of fluid and 
consequent extra refrigeration of tho blade as a whole. 

The advant3,go of methyl alcohol w'ill increase as the tem- 
perature is 1 o;>;er e d . and the effect of evaporation becomes 
less pronounced. 

So far , no a,ccount has been taken of the alcohol which 
is lost by evaporation. This can’ be calculated from equa- 
tion (S), but, since the coefficient of transfer of heat 
must be evaluated, the rate will be specific to a particu— 
1-ar position on the blade and particular velocity. Sy v/ay' 
of illustration, the- rate has been calculated for a blade of 
9— inch chord at the 5-per cent— chord station for -which ref- 
erence 5 gives the value of as 1.67. Eor a speed 

of 450 feet per second at .350 nilliraotors pressure, the val- 
ue -of kg, therefore, is 0.0016. The rate of eyaporation 

of ethyl alcohol from a correct alcohol— water mixture for a 
blade temperature of 23° E, from equation (3), is 0.53 pound 
per hour per sqxiare foot of area. The rate for pure ethyl 
alcohol, under the same conditions, is 3.13 poxinds per hour 
per square foot. 


The Carburetor 

The principal feature in the icing of cax-buretors of 
modern design is the defect in kinetic heating which occxirs 
when the velocity of the air is increased to above that of 
the free stream. This effect has already been discussed in 
connection with the temperature of the cambered surface of 
the blade of a propeller. In the case of the carburetor, it 
may be illustrated by a calculation of the tempera.ture of 
the surface w'hon the flo\>.- of air is throttled. 

Calculations have been for a teinporaturo at the entry 
to the carburetor of 38° E and a velocity of 938 feet por 
second. This temperature was chosen as being near to the 
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limit at which ice was observed to form in the tests de- 
scribed in reference 7 when no fuel was supplied to the car- 
buretor. The velocity has been calculated from Eernoulli's 
equation as that equivalent to the drop in: pressure across 
the throttle, given in table II of this reference — namely, 
10.7 inches of mercury. It is assumed that the air at in- 
let is saturated, that thei’e is no condensation, and that 
the flow is turbulent, Assuming the velocity at entry to, 
be zero , the value of t^, from equation (15) , is —34° E , 
and the value of Ati-, from equation (14), is 64.8° F. 

The air at entry is saturated with water vapor ; so, when 
throttled, the vapor pressure is 


—•ix X 5.8 = 3.64 millimeters 
29 


the barometric pressure at entry being assumed to be 29 
inches of mercury. If the v.'all is wotted v;ith water, equa- 
tion (10 ) gives the temperaturo at the surface t^ as 
28.7° F. Actually, at the wall thore will be. a gradient 
in temperature, so that heat conducted through’ the body of 
the carburetor v;ill cause the temperature to be somev/hat 
greater than that calculated. If ethyl alcohol is used to 
prevent freezing, the temperature at the surface f or . a cor- 
rect mixture , from equation (ll) , is 26.8° F, and 72.5 grams 
per 10 GO grams of water are rGqu.irGd.to depress the froe.z— 
ing point to' this temperature. 

At a.lt ittido , with the throttle fully open, the zone of 
high velocity is in the boost venturi in the Bendix— 
Stroraberg type of carburetor. Calculations have been made 
for this with air at entry to the carburetor at a tempera- 
ture of 34° F and pressure of 350 millimeters. A velocity 
of 600 foot per second was chosen, which gives a pressure at 
the throat of . 281 millimeters. The. velocity is higher, 
perhaps, than ordinarily would occur. On the assumption 
that no condensation occurs, the temperature at the, surface 
is calculated as 29.7° F, the method of calculation being 
the same as that discussed in the preceding paragraph. On 
the assumption that, condensation is complete, the tempera— 
ture at the surface is 32.9° F. 

These calculations are given by way of illustrating the 
method. Their primary use, it appears, is in calculating 
the rate of heating roquired to prot''ent ice by the method of 
references 1' and 2. The pfoblem in protection by fluids is 
to distribute the fluid efficiently. In this ca.se, the 

/ 
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method of calculation can he used to find the minimum rate 
of supply required, or to compare the merits of fluids of 
different composition, on the assumption of perfect distri- 
bution. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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TABLE I 

WET KIEETIC TEMPESATUHES EOH 600 ElET PER SECOITD, LAMIIAR PLOW 


Air temperature , 
(op) 

Pressure 
(mm ) 

Surface temperature , 

(op) 

A t In’ '? t 

A t d r 57 - 

0 

760 

19.4 

0.77 

0 

350 

CJl 

00 

.63 

25 

760 

40.1 

. 60 

35 

350 

36.5 

.49 


TABLE II .- TEMPERATURES AKD COWCIITTRAT lOlT OP ALCOHOL 
POR ALCOHOL-WATER MIXTURES OH PROPELLER BLADE 


Veloc- 
' ity 

(fps) 

Air 

tanper- 

ature, 

^0 

(°P) 

Pres- 

sure 

(mm) 

450 

15 o 2 

760 

450 

15.2 

760 

450 

19.2 

350 

450 

19.2 

350 

250 

24.2 

350 

250 

24.2 

350 

250 

19.2 

350 

250 

19.2 

350 

450 

15.2 

760 

450 . 

19.2 

350 

450 

1.6 

760 

450 

1.6 

760 

450 

5.5 

350 

450 

5.5 

350 

450 

21.9 

350 

450 

21.9 

350 

450 

19.2 

350 

450 

19*2 

350 


Su rf ace 
tempera- 
ture, to 


23.0 

26.3 


14.7 

11.5 


10,6 
15 , s 


19.9 

26 , s 


Cone on- • 


Excess 

t rat ion 


alcohol 

of fluid 
(grams/ 

Remarks 

required 

lOOOg 
water') • 


(percent) 

124 

92 

Por ethyl alcohol 
Por nonvolatile fluid 

'35 

124 

Por othjfl aleohol 

77 

70 

Por nonvolatile fluid 

124 

Por ethyl alcohol 

66 

75 

Por nonvolatile fluid 

1 S 4 

Por ethyl alcohol 

36 

135 

Por nonvolatile fluid 


Pure ethyl alcohol in 



air sat\iratGd with 
water vapot 


270 

Por ethyl alcohol 

13 

• 239 

Por nonvolatile fluid 

1 

270 

Por ethyl alcohol 

27 

212 

Por nonvolatile fluid 

S 6 

Por methyl alcohol 

21 s 

27 

Por nonvolatile fluid 

113 

Por methyl alcohol 

132 


Por nonvolatile fluid 

















